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Abstract The mean abundances of Mg, Si, Ca, Ti, Cr and Fe based on both strong and weak lines of 
Alpha Centauri A are determined by matching the observed line profiles with those synthesized from 
stellar atmospheric models and comparing these results with a similar analysis for the Sun. There is 
good agreement between the abundances from strong and weak lines. 

Strong lines should generally be an excellent indicator of abundance and far easier to measure than 
the weak lines normally used. Until the development of the Anstee, Barklem and O'Mara theory for 
collisional line broadening, the uncertainty in the value of the damping constant prevented strong lines 
being used for abundance determinations other than in close differential analyses. 

We found that Alpha Centauri A has a mean overabundance of 0.12±0.06 dex compared to solar mean 
abundances. This result agrees remarkably well with previous studies that did not use strong lines 
or the Anstee, Barklem and O'Mara theory for collisional line broadening. Our result support the 
conclusion that reliable abundances can be derived from strong lines provided this new theory for line 
broadening is used to calculate the van der Waal's damping. 

Keywords: stars: abundances, stars: individual (Alpha Centauri A), Sun: abundances 



1 Introduction 

The Alpha Centauri system has a special fascination 
for astrophysicists because it is the closest stellar sys- 
tem and its principle component, Alpha Centauri A (a 
Cen A), has a spectral type very similar to the Sun, 
G2V. a Cen A is also one of the brightest stars in the 
sky enabling spectra with a high spectral resolution 
and high signal to noise to be obtained. Most past 
analyses have concluded that the metal abundance of 
a Cen A is greater than that of the Sun. Analyses, 
such as Furenlid & Meylan (1990) covering 26 ele- 
ments from 500 lines and Neuforge-Verheecke & Ma- 
gain (1997) investigating 17 elements, concluded that 
the average metal overabundance compared to the Sun 
is 0.12±(0.02-0.04) dex and 0.24 dex respectively. An- 
other study, Chmielewski et al. (1992), concluded that 
a Cen A can be classed as a super metal rich star. 

The importance of the present study is the use of 
the Anstee, Barklem, and O'Mara 1 theory (ABO the- 
ory, Barklem et al. (1998a) ) to determine van der 
Waal's damping (VDW) for collisional line broadening 
for q Cen A. The ABO theory provides precise theo- 
retical damping constants (as demonstrated in recent 
results for the analysis of strong lines in the solar spec- 
trum by Allende Prieto et al.(2001)), which enables the 
use of strong lines for which reliable laboratory f- values 
exist. Strong line wings are also relatively insensitive 
to the effects of turbulence in the atmosphere. Thus 
strong line wings, together with the ABO theory, may 



1 Jim O'Mara was the primary supervisor for this 
project. Sadly he died suddenly in April 2002 during a 
trip to Italy and France before this work was finished. 



be used as reliable abundance indicators for elements 
where such strong lines exist. 

The choice of lines for this project follows Allende 
Prieto et al.(2001) investigation for the Sun and in- 
cludes weak neutral, weak ionised, and strong lines for 
six elements, Mg, Ca, Si, Ti, Cr, and Fe. 

We employ high quality coude echelle CCD spectra 
observed using the 74 inch telescope at Mt. Stromlo 
Observatory. The mass, distance, luminosity, and colours 
of a Cen A are used to determine its effective tem- 
perature, T e ff, and surface gravity, log(g s ) necessary 
to customise existing atmospheric models to match a 
Cen A. The model is based on Kurucz solar models, 
(Kurucz 1979). The customised model, along with rel- 
evant atomic data, are used to synthesize line profiles 
for a Cen A. Abundances are determined by match- 
ing the observed profiles, following determination of 
turbulence parameters. Steps are taken to verify the 
a Cen A model used. A solar model based on the 
Holweger-Miiller model (Holweger & Miiller 1974) is 
used to synthesize line profiles which are matched to 
the observational line profiles from the Jungfraujoch 
Atlas (Delbouille & Roland 1995) to determine the so- 
lar mean abundance. The solar and a Cen A mean 
abundances are compared to find the mean under- or 
overabundance for a Cen A. 



2 ABO Collisional Line Broad- 
ening Theory 

Collisional or VDW broadening is broadening resulting 
from the collision of atoms in the photosphere, ft is 
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especially important in cool stars such as our Sun and 
a Cen A that have predominately neutral hydrogen 
photospheres. This collisional broadening produces a 
Lorentzian line profile. The damping constant for col- 
lisional broadening is r co ;; and is included when syn- 
thesizing the line profiles. 

As discussed in Barklem et al. (1998a), up until 
the 1970's, formulations of VDW broadening devel- 
oped by Lindholm, Foley, and Unsold were used and 
it was widely held that a better theory was needed. 
These theories deal with VDW interactions between 
perturbing hydrogen atoms and the absorbing atom. 
Although the term van der Waal's broadening is still 
used today, it is really a misnomer as the actual line 
broadening theory is much more complex. 

K.A. Brueckner (Brueckner 1971) introduced a per- 
turbation theory formulation that involved long range 
interaction where the electron exchange could be ne- 
glected but not the overlap in the atomic charge distri- 
bution, a point taken up in O'Mara (1976). O'Mara's 
work deals with collision broadening theory and draws 
on work from many sources to develop the beginning of 
what has come to be known as the ABO theory of line 
broadening. This theory has been further developed 
by Anstee, Barkelm and O'Mara (Anstee & O'Mara 
1991, 1995; Barklem & O'Mara 1997; Barklem et al. 
1998b). Further development of the ABO theory and 
its relevances to solar and late-type star abundances is 
on-going with a code available on the world wide web 
(Barklem et al. 1998a) to calculate VDW. 

3 Observation and Data Re- 
duction 

Our spectra are taken on three separate observational 
runs in 1996 June/July and 2001 May on the 74 inch 
telescope at Mt Stromlo Observatory. The 120 inch fo- 
cal length coude camera is used with a 31.6 groove/mm 
echelle, cross dispersed with a 150 lines/mm grating. 
Several different wavelength settings of the echelle grat- 
ing and cross-disperser gratings are used to obtain al- 
most complete wavelength coverage from 4000-8000 A. 

The signal to noise ratio (S/N) varied with order 
across the CCD because of vignetting resulting from 
the cross-disperser not being near a pupil. The CCD 
has a gain of 2 electrons per ADU and most of the ex- 
posures are aimed at about 60000 electrons maximum. 
As the data are 60000-180000 electrons per resolution 
element, the nominal S/N ratio is 200-400. The ac- 
tual resolution was 125000 estimated by measuring the 
width (FWHM) of the line at wavelength 8252.379 A 
from the thorium arc spectrum: 



which corresponds to a velocity of 2.392 km s~ . 

The reduction process includes cleaning the raw 
spectra of cosmic rays and flat fielding. The spec- 
tral orders are extracted and any scattered light be- 
tween the orders is removed. Extracted spectra from 
a nearby "smooth spectrum" star /3 Centauri, are di- 
vided through the spectra of the target star, a Cen A 



to eliminate telluric (Earth's atmosphere) lines. An Th 
arc spectrum is used for wavelength calibration pur- 
poses. The spectra are smoothed to reduce noise and 
the continuum level is flattened. The wavelengths are 
also corrected for the radial velocity of a Cen A. 

4 The Atmospheric Models and 
Line Profiles 

Line profiles are synthesised for each line of interest 
and used to determine the abundances by fitting to 
the observed line profiles. Parameters such as macro 
and microturbulence, VDW damping, energy levels of 
the transitions, log gf values, and starting abundances 
are input to enable this direct matching process. 

Two models and subsequent line profiles are used 
for a Cen A. Model AK is an interpolated Kurucz 
model (Kurucz 1979) grid; model AH is a scaled Hol- 
weger & Miiller solar model (Holweger & Miiller 1974). 
Both models use the calculated values for a Cen A's 
effective temperature and surface gravity. These two 
models are compared to verify the validity of using the 
interpolated Kurucz model in this project. 

The third model SH is the Holweger & Miiller model 
(Holweger & Miiller 1974). This is used to synthesize 
line profiles to fit observed solar data obtained from 
the Jungfraujoch Atlas. 

The initial models based on published solar abun- 
dances, are used to compute line profiles to be com- 
pared with observations. If the line does not fit, the 
abundance of the element is adjusted and new number 
densities, opacities, and pressures are computed for a 
second iteration. This process is iterated to conver- 
gence yielding an abundance for each line. The tem- 
perature structure, Ti ogT , of the initial atmosphere is 
not adjusted. 

4.1 Input Parameters and Calcula- 
tions for the a Cen A Model At- 
mosphere and Synthesised Line 
Profiles 

4.1.1 Effective Temperature and Surface 
Gravity 

The parameters, effective temperature, T e //, and sur- 
face gravity, log(g s ), used to calculate the atmospheric 
model are very important. 

In Table we list the measurements and derived 
data for a Cen A and the Sun that are used to deter- 
mine the surface gravities and effective temperatures 
and are computed using the following equations: 

log.g, = log(^) (2) 

r ^=y^r (3) 

The derived values are in close agreement with 
those researched by other authors and are listed in Ta- 
blc|21 For a Cen A the resulting values are T e //=5784±5K 
and log(3 a )=4.28±0.01. 
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Table 1: Researched and Calculated* values used to customise the model atmosphere and produce synthesised line profiles 
for a Cen A's and the Sun 
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^'Blackwell & Shallis 1977; ( 2 'Perryman 1997; ( 3 >Ahrens 1995; ( 4 >Bessell et al. 1998; < 5 >Demarque et al. 1986, * calculated, 
error insignificant 



Table 2: Comparing parameters for a Cen A from reviewed and this paper 



Ref. 


T e //(K) 


log(gs) 


£(km s 1 ) 




L 


r 


m 


Alog(A)(dcx) 


1 


5784.3±5.5 


4.28±0.01 


1 


4.35±0.006 


1.57±0.3 


1.25±0.0024 


1.085±0.01 


0.12±0.06 


2 


5770±20 










1.23 






3 


5710±25 


4.0±0.2 


1.0±0.2 


4.38 




1.26 


1.085 


0.12±0.06 


4 


5765±50 








1.53 




1.085 


0.25±0.02 


5 


5800±20 


4.31±0.02 


1 


4.374±0.01 






1.085 




(i 


5830±30 


4.30±0.03 


1.09±0.11 




1.53 




1.085±0.01 





(^This paper; ( 2 >Soderblom 1986; < 3 )Furenlid & Meylan 1990; < 4 >Noels et al. 1991; ( 5 >Chmielewski et al. 1992; ( 6 )Neuforge- 
Verheecke & Magain 1997. 
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4.1.2 Input Data for Line Profile Synthesis 

The input data required for line synthesis are the en- 
ergy levels, log(gf), and VDW parameters taken from 
Allende Prieto et al. (2001) (Tables [3 & BJ , and the 
starting solar abundances from Grevesse & Sauval (1998) 
(Table 0. 

The parallax (tt) and angular diameter (8) are used 
to calculate a Cen A's radius: 



[d(0/36OO)(7r/18O)l 
r = j (5) 

The apparent magnitude is used to calculate the 
effective temperature with the following steps: 

M v = V-51og-^ (6) 
M bol = M V +BC (7) 



Microturbulence Values for the microturbulence 
are determined by using weak neutral, weak ionised, 
and strong Fe lines. The values fitted are 1.08±0.2 km 
s _1 for q Cen A with the AK model and 0.85 km s _1 
for the Sun using the SH model. 

4.3 Validity of the Model 

To judge whether the interpolated model AK is valid, 
comparison with an empirical model based on obser- 
vations is useful. 

The equivalent widths and profiles from model AK 
are compared to those produced with model AH a 
scaled solar empirical Holweger & Miiller model (Hol- 
weger & Miiller 1974). 

Eight lines are chosen and a good match of line 
profiles and equivalent widths are obtained. 

Another method for comparing the validity of model 
AK is to plot Ti ogT vs. logT for both models. As can 
be seen from Figure^there is close agreement for logr 
-3— 1-0.6. This is the range of optical depths that are 
covered in this research. A direct comparison can be 
seen in Figure |21 where a strong line, CrI at 4801A, for 
both models have been overlaid and match exactly. 
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4.2 Non-thermal Broadening Param- 
eters 

The non-thermal motions or turbulence that occur on 
a large scale compared with optical depth, is macro- 
turbulence, and on the small scale is microturbulence. 

Macroturbulence Macroturbulence will broaden 
the line profile but does not affect the line strength 
(equivalent width). The effects of stellar rotation (small) 
and instrumental profile are included as enhancements 
to the macroturbulence. 

Several single lines and one blend of lines are cho- 
sen to determine the effective macroturbulence. Var- 
ious values are tested on each single line to see if a 
common value can be used to match the general shape 
of the observed line. When the value for the effec- 
tive macroturbulence is determined, that value is used 
on the blended lines with no affect on the equivalent 
width of the synthesised line profile. 

The a Cen A value for the effective macroturbu- 
lence is determined to be 3.3 km s _1 . This includes 
2.4 km s" 1 due to instrumental effects and an atmo- 
spheric macroturbulence of 2.3 km s - . This indicates 
that most of the observed line broadening effects near 
the line core are from the large scale motions in the 
photosphere. Natural, Stark, and rotational broaden- 
ing are insignificant for the lines chosen for this study. 

The effective macroturbulence to match the solar 
model AH's line profiles to the observed solar data is 
1.6 km s" 1 . 



- Model AK 
■ Model AH 



AT vs.log T 



-3 -2.5 -2 -1.5 -1 -0.5 0.5 

log i 



Fi gure 1: Temperature vs. Optical Depth. Comparison 
of the AK and AH Models to test for for validity. 



5 Analysis 

5.1 Determining the Mean Abundance 

Fitting the shape of the synthesised line profiles from 
model AK against that of the observations is straight- 
forward for most of the lines. 

For all lines, care is taken to ensure that the area 
of the synthesised and observed line profiles are equiv- 
alent, allowing for differences in the exact matching of 
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Figure 2: Model AK & AH line profiles overlaid for CrI 
4801 A 



the line profile's core and wings. For two very weak 
lines and all the strong lines, extra care is needed. 

In the case of the strong lines, the core and wings 
cannot be matched simultaneously. In this case, the 
wings are matched, as per the ABO theory, and the 
core of the model's line profile is extended as low as 
possible. Care is taken that this technique is followed 
in both model's line profiles for those lines affected. 

As can be seen in Figure [3] the right hand wing of 
the line profile for the observed data for a Cen A is 
higher than that of the synthesised line profile. This 
does not happen in the Sun's line profile. The problem 
then is with the reduction process. Care is taken to 
ensure that features like this are compensated for in 
the determination of the abundances. 



5.2 Results 

The mean abundances are determined from 1 Mg, 7 
Ca, 9 Si, 13 Ti, 6 Cr, and 22 Fe lines. Each element is 
represented by species that cover weak neutral, weak 
ionised, and strong lines except for Mg that is repre- 
sented by one strong line, a Cen A shows a mean 
overabundance of 0.12±0.06 dex. The error is cal- 
culated from the variations in the individual species 
abundances. 

These abundances are larger than those used for 
starting abundances from Grevesse & Sauval (1998). 
The individual equivalent widths (W\), abundance (A), 
and A abundances (AA) for a Cen A and the Sun are 
listed in Tables El and H The mean abundance ({A)), 
the A mean abundances (A(j4,)) along with the start- 
ing abundances for a Cen A and the Sun for each el- 
ement and the final mean overabundance are listed in 
Tabled 



6 Summary and Discussions 

The mean abundance for the six elements investigated 
in the chemical composition of a Cen A are Mg=7.9 2 , 
Ca=6.46±0.09, Si=7.66±0.07, Ti=5.07±0.09, Cr=5.88±0.16, 
and Fe=7.71±0.20 dex. This leads to a mean abun- 
dance of 0.12±0.06 dex with respect to the Sun. 

Previous studies are not able to use strong lines as 
no reliable theory existed to calculate the collisional 
broadening of these lines. For this project the devel- 
opment of the ABO theory (Barklem et al. 1998a) to 
calculate the VDW damping, enable us to use strong 
lines in determining the mean abundance of a Cen A 
compared to the Sun. 

The a Cen A parameters used are T e //=5784.3±5.5K 
and log(g s )=4.28±0.01. Two models for a Cen A arc 
used, the second one (AH) a scaled solar model for 
comparison with the first (AK) to verify the validity 
of using an interpolated Kurucz model (Kurucz 1979). 
Once the use of model AK is validated, this model is 
used to synthesize line profiles for a Cen A to match 
the observed line profiles. 

Solar abundances are determined from comparing 
observed line profiles from the Jungfraujoch Atlas with 
those of the Holweger-Muller solar model (Holweger & 
Miiller 1974)). 

The results of this study, that a Cen A is over- 
abundant with respect to the Sun and can be included 
with other metal rich stars, agrees with those of pre- 
vious studies (Noels et al. 1991; Chmielewski et al. 
1992; Neuforge 1993; Neuforge-Verheecke & Magain 
1997) with an exact agreement of mean overabundance 
of 0.12±0.06 dex by Furenlid & Meylan (1990). This 
mean overabundance indicates that a Cen A did not 
originate in the same cloud as the Sun but from mate- 
rial that is more enriched by stellar nuclear processing. 

The above studies did not use strong lines or the 
ABO theory. Our result supports the determination 
that reliable abundances can be derived from strong 
lines provided the ABO theory is used to calculate the 
VDW damping. 

By using the ABO theory for strong lines, the anal- 
ysis of spectra, construction of model atmospheres, 
and the subsequent synthesised line profiles, the chemi- 
cal composition of apparently faint stars, such as those 
in external galaxies can be determined. All galaxies 
contain cool F,G K type stars whose spectra contain 
strong metallic lines. Previously these lines were not 
able to be used for absolute abundances due uncertain- 
ties in the theory for calculating the VDW damping. 
With the ABO theory now firmly established and re- 
liable model atmospheres existing for cool stars it will 
be possible to extend reliable abundance analyses to 
more and more distant galaxies using strong lines. 
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Table 3: Input Parameters and Abundances for a Ccn A and the Sun 
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Sil 


4.95 


-1.550 


1798 


0.221 


77.74 


7.610 


66.87 


7.500 


0.110 


5690.425 


Sil 


4.93 


-1.770 


1772 


0.222 


62.79 


7.616 


52.84 


7.520 


0.096 


5701.106 


Sil 


4.93 


-1.950 


1768 


0.222 


55.19 


7.680 


39.91 


7.510 


0.170 


5708.402 


Sil 


4.95 


-1.370 


1787 


0.222 


95.29 


7.640 


84.96 


7.537 


0.103 


5787.922 


CrI 


3.32 


-0.083 


1097 


0.291 


87.12 


6.203 


48.17 


5.660 


0.543 


5866.457 


Til 


1.07 


-0.784 


259 


0.262 


59.83 


5.120 


43.38 


4.980 


0.140 


5916.254 


Fel 


2.45 


-2.990 


341 


0.238 


64.02 


7.710 


52.80 


7.650 


0.060 


5922.115 


Til 


1.05 


-1.410 


313 


0.242 


28.10 


5.065 


17.82 


5.000 


0.065 


5948.545 


Sil 


5.08 


-1.130 


1875 


0.222 


107.51 


7.663 


94.74 


7.506 


0.157 


6082.715 


Fel 


2.22 


-3.570 


306 


0.271 


47.38 


7.697 


33.24 


7.590 


0.107 


6092.799 


Til 


1.89 


-1.323 


398 


0.239 


7.92 


5.125 


3.77 


4.96 


0.165 


6151.623 


Fel 


2.18 


-3.300 


277 


0.263 


62.73 


7.711 


49.26 


7.61 


0.101 


6161.297 


Cal 


2.52 


-1.266 


978 


0.257 


79.49 


6.540 


66.02 


6.409 


0.131 


6162.183 


Cal 


1.89 


-0.097 


878 


0.236 


280.96 


6.350 


277.82 


6.350 





6166.441 


Cal 


2.52 


-1.142 


976 


0.257 


81.87 


6.454 


75.14 


6.420 


0.034 


6173.342 


Fel 


2.2 


-2.880 


281 


0.266 


81.34 


7.740 


70.01 


7.664 


0.076 


6200.321 


Fel 


2.61 


-2.440 


350 


0.235 


85.81 


7.750 


74.75 


7.680 


0.070 


6258.109 


Til 


1.44 


-0.299 


355 


0.237 


62.19 


5.010 


51.87 


4.980 


0.030 


6297.801 


Fel 


2.2 


-2.750 


278 


0.264 


86.68 


7.700 


76.60 


7.650 


0.050 
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Figure 3: Fitting Model line profile to Observed Data for a Cen A and Sun to determine equivalent width and abundance 
for Strong Call line 8542.120 A (Scale: 1 divisions) 



Table 4: Input Parameters and Abundances for a Cen A and the Sun (Continued) 



A 








VDW< 2 > 


a Cen A' 1 ' 


Sunt 1 ) 






(A) 


Species^ 2 ' 


Elow (2) 


log(gf)(2) 


(7 


a 


W A (ml) 


A 


W A (ml) 


A 


(dcx) 


6371.361 


Sill 


8.12 


-0.000 


389 


0.189 


43.34 


7.850 


35.22 


7.520 


0.330 


6432.684 


Fell 


2.89 


-2.510 


174 


0.270 


53.76 


7.600 


43.42 


7.449 


0.151 


6455.604 


Cal 


2.52 


-1.290 


365 


0.241 


69.90 


6.490 


56.85 


6.350 


0.140 


6481.878 


Fcl 


2.28 


-2.980 


308 


0.243 


76.78 


7.757 


65.20 


7.680 


0.077 


6498.945 


Fcl 


0.96 


-4.700 


226 


0.253 


57.82 


7.740 


43.20 


7.658 


0.082 


6499.656 


Cal 


2.52 


-0.818 


364 


0.239 


101.06 


6.600 


92.17 


6.500 


0.100 


6516.086 


Fell 


2.89 


-3.380 


174 


0.270 


45.07 


7.270 


53.99 


7.550 


-0.280 


6750.161 


Fcl 


2.42 


-2.620 


335 


0.241 


111.89 


8.150 


79.07 


7.695 


0.455 


6978.861 


Fcl 


2.48 


-2.500 


337 


0.241 


108.08 


8.000 


82.04 


7.662 


0.338 


7357.735 


Til 


1.44 


-1.06 


329 


0.244 


29.57 


5.060 


20.77 


5.030 


0.030 


7515.836 


Fell 


3.9 


-3.450 


187 


0.271 


23.33 


7.750 


14.73 


7.517 


0.233 


7680.271 


Sil 


5.86 


-0.590 


2107 


0.495 


101.18 


7.640 


101.78 


7.610 


0.030 


7711.730 


Fell 


3.9 


-2.450 


186 


0.264 


58.42 


7.560 


50.79 


7.430 


0.130 


7918.387 


Sil 


5.95 


-0.510 


2934 


0.232 


108.85 


7.600 


102.00 


7.504 


0.096 


8327.067 


Fcl 


2.20 


-1.525 


258 


0.247 


226.21 


7.660 


198.38 


7.621 


0.039 


8542.120 


Call 


1.7 


-0.463 


291 


0.275 


3206.68 


6.360 


3274.47 


6.360 





8662.169 


Call 


1.69 


-0.723 


291 


0.275 


2568.95 


6.420 


2566.40 


6.380 


0.040 


8806.778 


Mgl 


4.33 


-0.120 


531 


0.292 


679.65 


7.900 


562.10 


7.680 


0.220 
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Tabic 5: Mean Abundances and Mean Overabundances for a Ccn A models (AK), and Sun Model (SH) 



Element Grevesse & Sauval 1998 a Cen A - (A) Solar - (A) A (A) (dex) 



Mg 


7.58±0.05 


7.9 (1 line) 


7.6800 


0.200 


Ca 


6.36±0.02 


6.459±0.092 


6.396±0.054 


0.064±0.017 


Si 


7.55±0.05 


7.664±0.0756 


7.521±0.0370 


0.144±0.011 


Ti 


5.02±0.06 


5.075±0.088 


4.983±0.032 


0.092±0.019 


Cr 


5.67±0.03 


5.880±0.163 


5.688±0.038 


0.192±0.029 


Fe 


7.50±0.05 


7.714±0.197 


7.599±0.098 


0.115±0.029 




a Cen A's mean 


overabundance = 


0.12±0.06 dex 





thank the anonymous referees for their constructive 
comments. 
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